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We report on our recent investigations in photonuclear production of heavy quarks and 
vector meson in ultraperipheral heavy ion collisions. In particular, our theoretical predic- 
tions are compared with the recent experimental measurements on coherent p (STAR) 
and J/*!> (PHENIX) photoproduction at RHIC and estimates for LHC are given. 
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1. Introduction 

In the last years, there has been a lot of interest in the description of electron- nucleus 
collisions at high energies. The results of current analysis show that future electron- 
nucleus colliders at HERA and RHIC, probably could determine whether parton 
distributions saturate and constraint the behavior of the nuclear gluon distribution 
in the full kinematical range 1,2 . However, until these colliders become reality we 
need to consider alternative searches in the current and/or scheduled accelerators 
which allow us to constraint the QCD dynamics. In this contribution, we summarize 
our investigations in Refs. 3 ' 4 about the possibility of using ultraperipheral heavy 
ion collisions (UPC's) as a photonuclear collider and study the heavy quark and 
vector meson production assuming distinct approaches for the QCD evolution. 

In heavy ion collisions the large number of photons coming from one of the col- 
liding nuclei will allow to study photoproduction, with energies W^a reaching to 
950 GeV for the LHC. The photonuclear cross sections are given by the convolu- 
tion between the photon flux from one of the nuclei and the cross section for the 
scattering photon- nuclei. The expression for the production of a given final state Y 
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(Y = QQ, V) ultraperipheral heavy ion collisions is then given by, 

oo 

(w 2 a = 2cj^S^) (1) 

where cu m in — My/4-fLirip, \/5nn is the c.m.s energy of the nucleus-nucleus system 
and jl is the Lorentz boost of a single beam. The Lorentz factor for LHC is 7l = 
2930, giving the maximum c.m.s. jA energy W 1 a < 950 GeV. The requirement 
that photoproduction is not accompanied by hadronic interaction (ultraperipheral 
collision) can be done by restricting the impact parameter b to be larger than twice 
the nuclear radius, Ra = 1.2 A 1 ^ fm. An analytic approximation for AA collisions 
can be obtained using as integration limit b > 2R A , producing 



dN (uu) 2Z 2 a em \ r? / , , . n , xN 

~itr = nJ n Ka {fj) Kl {fj) + 2 ^ {f]) ~ K ° {fj) ^ 



(2) 



where fj = IloRaHl and K 01 (x) are the modified Bessel functions. Therefore, the 
main ingredient for computing the production of the specific final state Y at UPC's 
is the information about its cross section in photon- nuclei interactions, g 1 a^yx , as 
defined in Eq. 1. In what follows, we perform phenomenology on heavy quarks and 
vector mesons production using the available high energy approaches. 

2. Photonuclear production of heavy quarks 

For our further analysis on photonuclear production of heavy quarks we will con- 
sider distinct high energy approaches, which are contrasted against the collinear 
approach 3 : (i) semihard formalism, (ii) saturation model (within the color dipolc 
formalism) and (iii) Color Glass Condensate (CGC) approach. The main goal is a 
comparison among those approaches and show whether the production at UPC's 
allows disentangle the underlying QCD dynamics at high energies. 

(i) Semihard formalism : In the fcj_-factorization (semihard) approach, the rel- 
evant QCD diagrams are considered with the virtualities and polarizations of the 
initial partons, carrying information on their transverse momenta. The scattering 
processes are described through the convolution of off-shell matrix elements with the 
unintegrated parton distribution, T(x, fcj_). Considering only the direct component 
of the photon, the cross section reads as (See, e.g. Ref. 5 ), 

„2 



semihard _ ' ' ' ' 



V-fA^QQX 



a. [ dz d* Pl± d* k± a ^™fo*L) a{z) ( at + ft-i-Pix) 

J fej_ ! V D\ D 2 




where D x = pf ± + rriQ and D 2 = (k± - p 1± ) 2 + m 2 Q and A(z) = [z 2 + (1 — 
z) 2 ]. The transverse momenta of the heavy quark (antiquark) are denoted by p 1± 
and p 2± = (fej_ — Pu_), respectively. The heavy quark longitudinal momentum 
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fraction is labeled by z. For the scale fx in the strong coupling constant we use the 
prescription [i 2 = k\ + nig. We use also the simple ansatz for the unintegrated 

gluon distributions, — dxGA ( x _\k±) where xGa(x, Q 2 ) is the nuclear gluon 

9 In K ± 

distribution (see Ref. 3 for details in the numerical calculations). 

(ii) Saturation model : Based on the color dipole formalism, it is an extension 
of the ep saturation model through Glauber-Gribov formalism. In this model the 
cross section for the heavy quark photoproduction on nuclei targets is given by 6 ' 7 , 

<TW = [ dz J d * r \*t(z, r, Q 2 = 0)| 2 r 2 ,A), (4) 

where the transverse wave function is known 7,8 . The nuclear dipole cross section is 
given by 6 ' 7 , 



<4,(£, r 2 ,A)= J d% a|l - ex p^~AT A (b)a p dip (x, r 2 )] J 



(5) 



where b is the impact parameter of the center of the dipole relative to the center of 
the nucleus and the integrand gives the total dipole-nucleus cross section for fixed 
impact parameter. The nuclear profile function is labeled by Ta(V). The param- 
eterization for the dipole cross section takes the eikonal-likc form, <J dip (x, r 2 ) = 
do [1 — exp (— Q 2 {x) r 2 /4) ], where one has used the parameters from saturation 
model, which include the charm quark with mass m c = 1.5 GeV and the definition 
x = (Q 2 +4m 2 ? )/W 2 j4 . The saturation scale Q 2 s (x) = {x /x) X GeV 2 , gives the 
onset of the saturation phenomenon to the process. The equation above sums up 
all the multiple elastic rescattering diagrams of the qq pair and is justified for large 
coherence length, where the transverse separation r of partons in the multiparton 
Fock state of the photon becomes as good a conserved quantity as the angular mo- 
mentum, namely the size of the pair r becomes eigenvalue of the scattering matrix. 

(iii) Color Glass Condensate : The regime of a CGC is characterized by the 
limitation on the maximum phase-space parton density that can be reached in the 
hadron/nuclear wavefunction (parton saturation) and very high values of the QCD 
field strength F^ v w l/a s (For a review see, e.g, 1 ). The large values of the gluon 
distribution at saturation suggests the use of semi-classical methods, which allow 
to describe the small-x gluons inside a fast moving nucleus by a classical color 
field. In Refs. 9 the heavy quark production in UPC's has been analyzed in the 
CGC formalism. In Ref. 3 , we have improved that analysis using a realistic photon 
flux and a color field correlator including quantum radiation effects. The input 
photoproduction cross section reads as 3 , 



4(fcj - m 2 Q ) 



where we define the rapidity Y = ln(l/x) = ln(2w7i/4mg) and e q is the quark 
charge. In Ref. 3 , we obtained the following analytical expression for the color field 
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Table 1. The photonuclear heavy quark cross sections for UPC's at LHC. 



QQ 


Collinear 


SAT-MOD 


SEMIHARD I (II) 


CGC 


cc 


2056 mb 


862 mb 


2079 (1679.3) mb 


633 mb 


bb 


20.1 mb 


10.75 mb 


18 (15.5) mb 


8.9 mb 



correlator considering a suitable ansatz for the dipole-nucleus cross section, 
~ 47r / k 2 \ 

C{x > k±)= m^ exp \-QM)' (6) 

where we have assumed for the nuclear saturation scale, Q 2 sA {x) — A 1 / 3 x Q 2 s {x). 

Having summarized the theoretical models for heavy quark photoproduction, 
let us present the numerical calculation of their total cross section at UPC's. We 
focus mostly on LHC domain where small values of x would be probed. In the 
following, one considers the charm and bottom masses m c = 1.5 GeV and mi, — 
4.5 GeV , respectively. Moreover, for PbPb collisions at LHC, one has the c.m.s. 
energy of the ion-ion system V<Snn = 5500 GeV. The results are presented in Table 
1. The collinear approach gives a larger rate, followed by the semihard approach 
(labels I and II refer to GRV94 and GRV98 gluon pdf's). The saturation model and 
CGC formalisms give similar results, including a closer ratio for charm to bottom 
production. Concerning the CGC approach, our phcnomcnological educated guess 
for the color field correlator seems to produce quite reliable estimates. Therefore, 
the photonuclear production of heavy quarks allow us to constraint already in the 
current nuclear accelerators the QCD dynamics since the main features from photon- 
nuclei collisions hold in the UPC reactions. 



3. Photonuclear production of vector mesons 

Let us consider the scattering process -fA — > V A in the QCD dipole approach, 
where V stands for both light and heavy mesons. The main motivation for using 
this approach is the easy intuitive interpretation and its successful data description 
in the proton case. The scattering process can be seen in the target rest frame as a 
succession in time of three factorizable subprocesses: i) the photon fluctuates in a 
quark-antiquark pair (the dipole) , ii) this color dipole interacts with the target and, 
iii) the pair converts into vector meson final state. In the color dipole formalism, the 
corresponding imaginary part of the amplitude at zero momentum transfer reads 
as 12 < 13 , 



XmA 1 A^vA = 




where ^ ^(z, r) and ^^(z, r) are the light-cone wavefunctions of the photon and 
vector meson, respectively. The quark and antiquark helicities are labeled by h and 
h and reference to the meson and photon helicities are implicitly understood. The 
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A STAR DATA 
CUTS: t<0.02 

lyl<3 



APHENIX(preliminarv|t 
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(b) 



Fig. 1. (a) Energy dependence of coherent p photoproduction at RHIC. (b) Rapidity distribution 
of coherent J/NE" photoproduction (with nuclear breakup) at RHIC (see text). 



variable r defines the relative transverse separation of the pair (dipole) and z (1 — z) 
is the longitudinal momentum fractions of the quark (antiquark). 

In the dipole formalism, the light-cone wavefunctions ^ h n (z, r) in the mixed 
representation (z,r) can be completely determined using light cone perturbation 
theory. On the other hand, for vector mesons, the light-cone wavefunctions are not 
known in a systematic way and they are thus obtained through models (For a re- 
cent detailed discussion see Ref. 11 ). Here, we follows the analytically simple DGKP 
approach 14 , which assumes that the dependencies on r and z of the wavefunction 
are factorised, with a Gaussian dependence on r. We keep the original parameters 
of the model. The main shortcoming of this approach is that it breaks the rota- 
tional invariance between transverse and longitudinally polarized vector mesons 11 . 
However, as it describes reasonably the HERA data for vector meson production, 
as pointed out in Ref. 15 , we will use it here. 

The corresponding parameters for the meson wavefunctions are presented in 
Table 1 of Ref. 10 . Following Ref. 13 , we have estimated contribution from real part 
for the photoproduction of vector mesons: it is about 3% for light mesons and it 
reaches 13% for J 10 . Additionally for heavy mesons we have take into account 
the skewedness effects, associated to off-forward features of the process, which are 
increasingly important in this case. Finally, the photonuclear cross section is given 

by 



(T 



tot 

jA^VA 



[2m A-yA- 



■VA 



16tt 



(1 + /3 2 ) 



dt\F(t)f 



with t m i n = (niy/2w) 2 . The quantity fi is the ratio between the imaginary and real 
part of the amplitude. We have used an analytical approximation 16 of the Woods- 
Saxon distribution as a hard sphere, with radius Ra, convoluted with a Yukawa 
potential with range a = 0.7 fm in order to compute the nuclear form factor, F(t). 
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Table 2. Cross sections for coherent meson photoproduction at LHC energy. 





HEAVY ION 


J/* (3097) 


0(1019) 


U! (782) 


p(770) 


LHC 


CaCa 


436 ub 


12 mb 


14 mb 


128 mb 




PbPb 


41.5 mb 


998 mb 


1131 mb 


10069 mb 



Recently, the STAR Collaboration at RHIC published the first experimental 
measurement of coherent p production in gold-gold UPC's at yfs = 130 GeV 
17 . The energy dependence of the cross section is presented in Fig 1-a. Our the- 
oretical prediction in the curve take into account the experimental cuts, which 
gives cr t hoory(|j/| < 3) = 410 mb, in good agreement with the STAR measurement 
0-star(|2/| < 3) = 370 ± 170 (stat) ± 80 (syst) mb. Furthermore, the PHENIX Col- 
laboration has preliminary measurements of the cross section for coherent J '/\& 
photoproduction in UPC at midrapidity at \/s = 200 GeV 18 , accompanied by nu- 
clear breakup. The rapidity distribution is shown in Fig 1-b in comparison with our 
prediction. The theoretical estimation gives da / dy\ y= o — 58 pb, which is in agree- 
ment with the PHENIX measurement dapHENix/dy\ y= o = 48 ± 16 (stat) ± 18 (syst) 
pb. Finally, in Table 2 one shows the predictions for the integrated cross sections for 
LHC. The experimental feasibility and signal separation on the reaction channels 
presented here are reasonably clear, namely applying a (low) transverse momentum 
cut pr < 1 and two rapidity gaps in the final state for the meson case. In contrast, 
for heavy quarks we have only one rapidity gap. 
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